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We discuss femtosecond Raman type techniques to control molecular vibra- 
te tions, which can be implemented for internal state cooling from Feshbach states 
^ ■ with the use of optical frequency combs with and without modulation. The 
'. technique makes use of multiple two-photon resonances induced by optical fre- 



quencies present in the comb. It provides us with a useful tool to study the 
details of molecular dynamics at ultracold temperatures. In our theoretical 
model we take into account decoherence in the form of spontaneous emission 
and collisional dephasing in order to ascertain an accurate model of the pop- 
ulation transfer in the three-level system. We analyze the effects of odd and 
even chirps of the optical frequency comb in the form of sine and cosine func- 
tions on the population transfer. We compare the effects of these chirps to 
the results attained with the standard optical frequency comb to see if they 
00 1 increase the population transfer to the final deeply bound state in the pres- 

o: 

CN ■ ence of decoherence. We also analyze the inherent phase relation that takes 

place owing to collisional dephasing between molecules in each of the states. 
^ ■ This ability to control the rovibrational states of a molecule with an optical 

frequency comb enables us to create a deeply bound ultracold polar molecules 
from the Feshbach state. 
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Decoherence at ultracold temperatures is the subject of particular interest and impor- 
tance in light of the development of the methods to manipulate ultracold gases, create and 



control ultracold molecules |2( , study ultracold collisions and chemical reactions [3| . De- 
coherence is inherently present in ultracold dynamics and we study it semiclassically within 
the process of creation of diatomic KRb molecules from Feshbach, weakly bound states. 
Experimentally ultracold polar KRb molecules were obtained using stimulated Raman adi- 
abatic passage (STIRAP), [4]. As a viable substitute to the STIRAP process, we created a 
method which makes use of optical frequency combs to perform two-photon resonances and 
transfer population in a stepwise manner from the Feshbach state to the ground, ultracold 
state. The implementation of an optical frequency comb (OFC) may be beneficial owing 
to its intrinsic ability to address the manifold of excited states simultaneously. We make 
use of a standard OFC, and the one with the sinusoidal phase modulation to induce two- 
photon Raman transitions. Our theory demonstrates that the impact of decoherence may 
be minimized by implementing the sinusoidally modulated OFC. 

The standard OFC is generated by the phase locked pulse train of the form 

E(t) = Sf- 1 E exp(-(t - kTf/(2T 2 )) cos (cu L t - kuj L T + 0). (1) 

Eq.flTJ gives a periodic envelope of the field which oscillates with an optical carrier frequency 
ul, and k is an integer number. The pulse train period T is much greater than r, where r 
is an individual pulse duration. A strictly periodic envelope function can be expressed as 
a spectrum by the Fourier series Y^^L^Aqexpi^— i(u>L + qco r )t) + c.c, where co r = 2tcT~ 1 is 
the pulse repetition rate and A q are spectral intensities that do not depend on time. The 
spectrum is a comb of laser frequencies precisely spaced by the pulse repetition rate u r . A 
qualitative picture of the temporal field and the frequency comb are shown on Fig.((T]a,b). 
The period of the pulse train may be within a 10 ns time scale, e.g. js], it determines the 
spacing between modes in the frequency comb to be 100 MHz. The Ti:Saphfire laser gives 
rise to about 10 6 modes in the comb. 

We study the phase variation across an individual pulse in the form of the sin function 

E(t,z) = l/2^~ 1 E exp(-(t-kT)y(2T 2 ))cos(uj L (t-kT) + ^osm(n(t-kT)) + (j)). (2) 

Here the $ is the modulation amplitude and the Q is the modulation frequency. The time- 
dependent phase of the field in the form of a sin function brings additional peaks to the 




Fig. 1. (a) A qualitative figure of the pulse train in Eq.(p]and (b) the respective frequency 
comb, (c) The sinusoidally modulated pulse train as in Eq.(J2]) and (d) the frequency comb 
with sets' uniform spacing determined by Q and spacing between the modes precisely equal 
to 2-kT- 1 . 

frequency comb. Laser frequency ujl determines the center of the frequency comb, while Q 
determines the uniform spacing of the individual sets of modes. Within each set of modes, 
the uniform spacing between the modes is determined by 2irT~ l . Eq.(J2D gives the pulse train 
and the intensity spectrum qualitatively shown in Fig.fll] c,d). Increasing Q broadens the 
frequency comb region, while $o determines the amplitude of the comb's peaks. Frequency 
comb peaks are narrower for longer propagation time of the pulse train, that is l/t to tai 
determines the bandwidth of an individual mode. 

A sinusoidal modulation at the MHz frequency of the beam of a cw ring dye laser was 
implemented in [6] to study absorption resonances in I 2 with high precision. Here we assume 
iz region for which the approach which makes use of an enhancement cavity 

We consider a three-level A system interacting with a frequency comb do describe semi- 
classically the dynamics of the Feshbach molecules interacting with radiation resulting in 
their stepwise transformation into the ultracold molecules. The A system is formed by the 
energy level that include the Feshbach state, the excited, transitional state and the final, 
ultracold state, Fig(J2J). 

For the description of the system dynamics we refer to the Leouville-von Neumann equa- 



modulation in T 
may be applied 
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Fig. 2. A three-level A system modeling the Feshbach state, the excited state manifold and 
the final state of molecules interacting with a phase locked pulse 

tion which provides a set of coupled differential equations for density matrix elements 

p u = 2Im[H 12 p2i + H 13 p 31 ] 

p 22 = 2Im[H 2 iPu + H23P32} 
p 33 = 2Im[H 31 p 13 + H 32 p 23 ) 

P12 = -iH 12 (p22 - Pu) - iHi 3 p 32 + iH 32 pi3 (3) 
pi 3 = -iH 13 (p 33 - pu) - iH 12 p23 + iH 2 3Pi2 
P23 = ~iH 23 {p 33 - P22) - iH 2 iPi 3 + iH 13 p 2 i 

The interaction Hamiltonian used in these equations is written beyond the rotating wave 
approximation. Nonzero Hamiltonian matrix elements in the interaction representation read 
H ji = n R (t-T)[exp{-i((u L +u ji )(t-T)+M(t-T))}+ex^ 

here, i,j are the indexes of the basis set, i — 1, 2 and j + M(t — T) = $0 sin Q(t — T) is 
the phase modulation in a single pulse which is zero for the standard OFC, and Quit — T) = 
fi^exp (— {t — T) 2 /(2r 2 )) is the Rabi frequency, with the peak value Qr. The decoherence 
effects caused by spontaneous emission and collisions are taken into account through the 
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reduced density matrix elements 

Pn)s P = 721P22 Pi2) sp ,coi = -{^y + 1 ¥ + r 2 i)pi2 

f>22)sp = -721P22 - 723P22 Pn)sp,col = -r 3 lPl3 ( 4 ) 

P33) S p = 723P22 P23)sp,cal = -(^f- + J f + r 2 3)P23, 

here 721,732 are spontaneous emission rates from the excited state |2 > to the initial state 
|1 > and the final state |3 > respectively. The 7s may also be attributed to inelastic collision 
rates. The collisional dephasing is taken into account by the rates r 2 i, r 3 i, r 23 . The reduced 
density matrix elements were added to Eqs. ([3]) which were solved numerically. 

To demonstrate the powerful tool of a modulated OFC to control rovibrational degrees 
of freedom at ultracold temperatures in the presence of decoherence, we analyze the impact 
of different channels individually within our model. At first, we discuss the case when deco- 
herence is not taken into account and emphasize similarities and differences in the dynamics 
induced by the standard and the modulated OFCs. Then, we turn on the decoherence in our 
model experiment, and illuminate the advantages the sinusoidally modulated OFC brings 
to ultrafast dynamics by minimizing the effects caused by fast spontaneous emission and 
collisions. 

The standard OFC and the sinusoidally modulated OFC application to a three-level A 
system leads to a full population transfer to the ultracold state from the Feshbach state s| . 
The difference is in the degree of involvement of the intermediate state into dynamics. The 
standard OFC populates substantially the excited, transitional state; its population may 
rise up to 50% depending on whether the carrier frequency of the pulse train is in resonance 
with the transition from the excited to ultracold state or is detuned. The results for the 
case of the one-photon resonance condition are presented on Fig.Q. 

In contrast to the effect induced by the standard OFC, the sinusoidally modulated OFC, 
while performing a stepwise population transfer to the final, ultracold state, only negligibly 
populates the excited state. This dynamics is observed under the condition of the one-photon 
resonances of the carrier and the modulation frequencies with the Feshbach-to-excited state 
and the excited-to-ultracold state transitional frequencies respectively. The dynamics is 
presented in Fig. (HI). 

Before discussing the impact of decoherence, let us point our attention to the dependence 
of the dynamics on the parity of the applied chirp of the field. The implementation of 
the odd and even chirps may result in a substantially different response of the system and 
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Fig. 3. Population dynamics in the three-level A system using an optical frequency comb 
having f r =5 GHz, and zero offset frequency. Parameters of the pulse train are the carrier 
frequency o;l=434.8 THz, the pulse duration r =3 fs, the peak Rabi frequency 0^=1.26 
THz; the system parameters are W2i=309.3 THz, W32=434.8 THz, and W3i=125.5 THz (J). 
Black line shows population of the ground state, red line - population of the transitional 
state, and green line - that of the final state. Time is given in the units of [w^ 1 ], where 
uj = (X>3i=125.5 THz. 

the quantum yield. We investigated the effect of the sine and cosine modulation within 
the three-level model and revealed an significant dependence of the population dynamics 
on the parity of the phase modulation. A sinusoidal modulation across an individual pulse 
in a phase locked pulse train leads to a stepwise adiabatic population transfer as described 
above, while a cosine modulation induces population jumps from the initial to the final state 
within a much shorter timescale and is accompanied by a significant population of the excited 
state during the transitional time. The phenomenon strongly depends on the magnitude of 
the Rabi frequency. According to [9^ the different results of the sine and cosine chirp are 
expected. The authors analyzed the Taylor series expansion of the field phase having an 
arbitrary time dependence and showed that even chirps lead to self-induced transparency, 
whereas odd chirps to population inversion. 

In the presence of fast decoherence, different degree of the involvement of the excited 
state into population dynamics induced by the standard and the sinusoidally modulated 
OFCs leads to a remarkably different response of the three-level system to the applied field. 
We decompose the contributions from spontaneous decay and collisional dephasing in our 



6 



10 20 30 

10 4 tca 



Fig. 4. Population transfer in the three-level A system, achieved via the resonant Raman 
transitions using a phase- modulated optical frequency comb described by Eq.([2j). The values 
of the system parameters are W32 = 410.7 THz, U21 = 340.7 THz, [10]. The carrier frequency 
ul = U32, the modulation frequency Q = u; 2 i, the modulation amplitude $0 = 4 and the peak 
Rabi frequency Qr=70 THz. Stepwise, adiabatic accumulation of the population is observed 
in state |3 >, (green), which is the ultracold KRb state. The population of the Feshbach 
state |1 >, (black), comes gradually to zero, while the excited state manifold |2 >, (red), is 
slightly populated during the transitional time. Full population transfer is accomplished in 
109 pulses. Time is given in the units of [w" 1 ], where u = U3i=70 THz. 

analysis. The time-dependent dynamics induced by the standard OFC in the presence of 
only collisional dephasing is presented by solid lines in Fig.flS]). Here, the quantum yield 
to the ultracold state is only 38% in the presence of fast collisions. When we add the 
spontaneous decay from the excited state, we observe that it reduces the collisional impact 
and yields to a higher steady state population of the final, ultracold state equal to 45%, see 
Fig.flSJ), (dashed lines). The data are presented for the following decoherence parameters: 
solid lines for 721 = 0., r 12 = T 2 3 = 0.001; and dashed lines for 721 = 732 = 0.001, and 



T 12 = T 2 3 = 0.001, in the units of u 3 i = 125THz, |4j. The Rabi frequency £Ir is equal to 0.1. 
The interplay of the spontaneous decay rate and the collision dephasing rate may result in 
an increase of the quantum yield of ultracold molecules. 

Next, we demonstrate the effect caused by the sinusoidal modulation across an individual 
pulse in the phase-locked pulse train. It results in a remarkably high quantum yield to the 
final, ultracold state under the same decoherence conditions. A comparison is shown on 
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Fig. 5. Population dynamics induced by standard optical frequency comb in the presence 
of spontaneous decay and collisions. Parameters of the pulse train are the carrier frequency 
wl=434.8 THz, the pulse duration r =3 fs, the peak Rabi frequency f2#=12.6 THz; the 
system parameters are o>2i=309.3 THz, ^32=434. 8 THz, and ^31 = 125.5 THz. Solid lines 
represent the case of pure collisional dephasing, r 31 = 0, r 2 i = r 32 = 0.001, dashed lines 
correspond to the same values of T's and 721 = 732 = 0.001. 

Fig. ([6]), where the dashed lines correspond to the standard OFC implementation same as 
in the previous figure, solid lines show the system response to the sine modulated OFC. The 
sinusoidally modulated comb negligibly involves the excited state into transitional dynamics 
(Fig.(J4])), minimizing the effects caused by decoherence. Two horizontal solid lines coinciding 
at the population equal to 0.5 show the steady state solution when the cosine modulation 
across an individual pulse in the pulse train was implemented justifying the importance of 
the parity of the chirp in achieving the desired quantum yield. Notably, we have uncovered 
a rigor expression for the relation between collisional dephasing rates in the three-level A 
system. Previously it was observed that for certain values of collisional dephasing rates the 
received populations at times incorrectly exceeded the initial total population. However, 
if to look from the physical point of view at several special cases, e.g., when molecules in 
ultracold configuration collide elastically with the atomic Rb and K at the same rate as the 
Feshbach molecules, then one would see that the state |3 > and |1 > amplitudes loose their 
phase synchronously and thus the r 31 must be zero, while T 2 i and r 32 must be not and 
be equal to each other. This example demonstrates that the values of collisional dephasing 
rates are mutually dependent. With the help of the phase diagram, we obtained a general 
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relation between dephasing rates in the three-level A system which reads 1^3 = Ti2 + I^i. 
With the values of collisional dephasing rates limited by this relation we observed that 
our total population always remained constant. Our results complement with the work in 



In summary, we have analyzed the interaction of an optical frequency comb with a three- 
levle A system that leads to cooling of internal degrees of freedom in the KRb molecule from 
the Feshbach state into ultracold deeply bound state. Using the reduced density matrix to 
introduce decoherence in the forms of spontaneous emission and collisional dephasing, we 
studied the effectiveness of the population transfer into the ultracold state performed by the 
standard OFC and a sine or cosine- chirped OFC and uncovered a remarkable difference in 
the system's response. We have shown that in the presence of fast decoherence an efficient 
cooling of internal degrees of freedom form the Feshbach state may be accomplished by a 
sinusoidally modulated OFC. While a standard OFC implementation leads to a substantial 
loss of the efficiency in the presence of fast dephasing. The sine modulated OFC induces 
the response in the three-level system that resembles the dark state formed in STIRAP. It 
negligibly populates the excited state manifold during a stepwise population transfer from 
the initial to the final, ultracold state and, thus, minimizes the effect of decay and optimally 
resolves the collisional impact. 
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